Quantification of memory effects in the spin-boson model 
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Employing a recently proposed measure for quantum non-Markovianity, we carry out a systematic 
study of the size of memory effects in the spin-boson model for a large region of temperature and 
frequency cutoff parameters. The dynamics of the open system is described utilizing a second-order 
time-convolutionless master equation without the Markov or rotating wave approximations. While 
the dynamics is found to be strongly non-Markovian for low temperatures and cutoffs, in general, 
we observe a special regime favoring Markovian behavior. This effect is explained as resulting 
from a resonance between the system's transition frequency and the frequencies of the dominant 
environmental modes. We further demonstrate that the corresponding Redfield equation is capable 
of reproducing the characteristic features of the non-Markovian quantum behavior of the model. 



I. INTRODUCTION 

Almost any realistic physical system interacts with its 
environment. The effects of this interaction can be de- 
scribed within the framework of open quantum systems 
[TJ The coupling of a quantum system to an envi- 
ronment typically leads to dissipation and decoherence 
processes in the open system. However, it is known that 
due to the interaction, non-Markovian effects can emerge 
as well. Following the classical definition [31 0], quantum 
non-Markovianity is often loosely described as the occur- 
rence of memory effects in which the environment acts as 
a memory allowing the earlier open system states to have 
an effect on the later dynamics of the system. It is known 
that non-Markovian effects can play an important role in 
many applications 3, 5] [6], but the absence of a proper 
definition of non-Markovianity in the quantum case im- 
peded systematic studies of these effects. Recently, sev- 
eral formal definitions for quantum non-Markovianity [Ti- 
ll], which are based on different mathematical and phys- 
ical concepts, have been proposed. 

Here, we employ the definition developed in Ref. [7J, 
which is based on the exchange of information between 
the system and its environment and enables a unique 
quantification of the size of quantum memory effects in 
the dynamics of a quantum system. We will study the 
behavior of this measure for a selected example: the 
spin-boson model, which is the paradigm of a dissipa- 
tive two-level system [21 [10] . The spin-boson model has 
been studied extensively over the past decades and has 
been successfully used to describe, e.g., chemical reac- 
tions, especially charge transfer processes in the con- 
densed phase [11], tunneling of defects and impurities 
in solids and metals [12) . and tunneling in amorphous 
materials [2] . Also, the famous Kondo problem [TS] 
could be linked to the spin-boson model [H], which led 
to the explanation for the Kondo effect [TTl [18] . Further, 
the dynamics of magnetic flux trapped in a SQUID ring 



has also been understood with the help of this model 
and, based on this setup, quantum computing devices 
have been proposed and partially implemented [T5H22]- 
The pursuit of the experimental realization of stable and 
controllable qubits has become highly important, and it 
is thus not surprising that quantum dots made of semi- 
conductors have also been investigated as possible qubit 
implementations [231 HI] ■ Such quantum dots and their 
dissipative interaction with the surrounding components 
are often described by the spin-boson model. Another 
very promising approach to quantum computation lies 
in trapped ions [25] . as these systems allow for a very 
high degree of control and long qubit lifetimes [23 [2Z] • 
Again, the spin-boson model is a key tool to a theoretical 
understanding of the important effects [28] . 

It is not only interesting to characterize memory ef- 
fects in this prototypical example, but an important mo- 
tivation for the present study also lies in the need for a 
better understanding and characterization of the prop- 
erties of the non-Markovianity measure itself, especially 
for more complex physical systems. There are several 
publications studying this quantity in more detail |29j . 
applying it to physical models [30U33] . and comparing it 
to other measures [3IH3E]- Additionally, several experi- 
mental measurements of the non-Markovianity measure 
have been carried out recently [33 [31], demonstrating 
that it provides an experimentally accessible observable 
which quantifies memory effects. Quite recently further 
theoretical implications and experimental applications to 
nonlocal quantum memory effects have been developed 
[39]. 

The paper is organized as follows: In Sec.[H]we briefly 
present the non-Markovianity measure of Ref. [7j , intro- 
duce the spin-boson model and the corresponding second- 
order time-convolutionless master equation, and explain 
the numerical simulation method used to determine the 
measure for quantum non-Markovianity. The numerical 
simulation results are discussed and linked to a theoret- 
ical interpretation in Sec. 



Ill where we also revisit the 
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role of the standard Markov approximation (which is not 
used in Sec. [Tl]). Finally, in Sec. IV 
results and draw some conclusions. 



we summarize our 
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II. QUANTIFYING NON-MARKOVIANITY 

A. Measure for the degree of memory effects 

We quantify memory effects in the dynamics of the 
spin-boson system, employing a recently proposed gen- 
eral measure for the degree of non-Markovian behavior 
in the evolution of an open quantum system [7|. This 
measure is based on the idea that Markovian time evo- 
lutions are characterized by a continuous loss of infor- 
mation from the open system to the environment, while 
non-Markovian dynamics feature a flow of information 
from the environment back to the open system. An ap- 
propriate tool to measure this information flow is given 
by the trace distance between two quantum states p\ and 
P2, which is defined by [40] 



2 Tr|Pl 



P7\ 



(1) 



where the modulus of an operator A is given by \A\ = 
V A. The trace distance can be interpreted as the dis- 
tinguishability of the states p\ and p 2 55 ■ We suppose 
that the dynamics of the open system is described by a 
family of trace preserving and completely positive maps 
Any pair of initial states pi,2(0) then evolves into 
Pi,2{t) = < &(t)pi,2(0) at time t > 0. If the family of maps 
$(t) is divisible, e.g., if it forms a semigroup with a gen- 
erator in Lindblad form 021 |43] , the trace distance 



V{t) =V( Pl (t),p 2 (t)) =X>(*(f) Pl (0), *(<)/*(<))) (2) 

between the time-dependent pair of states decreases 
monotonically [23]. This implies that the states p\ and 
P2 tend to become less and less distinguishable over time, 
which means that information about the system states is 
lost to the environment. On the other hand, whenever 
the trace distance increases, the two states become more 
distinguishable, and hence, information must have flowed 
back from the environment to the system, which is a clear 
signature for the presence of memory effects. 

On the basis of this interpretation, the measure for the 
non-Markovianity of the dynamics of an open system is 
defined by 



= max / dta(t), 

Pl,2(0) J (T> o 



(3) 



where a(t) — -^T>(t) denotes the rate of change of the 
trace distance given by Eq. ([2| and the time integra- 
tion is extended over all intervals in which a is positive. 
The measure, thus, quantifies the total backflow of in- 
formation from the environment to the open system. It 
further involves a maximization over the initial pair of 
states pi,2(0) and therefore represents a functional of the 
family of dynamical maps describing the physical 
process. 

By definition, the measure A/"($) is non-negative, and 
we have Af($) — if and only if the process is Marko- 
vian. A nonzero measure, J\f(<fr) > 0, implies that the 



process <!>(i) is nondivisible and cannot be described, for 
example, by a time-local master equation with positive 
rates [22] ■ We note that the non-Markovianity measure 
AT(Q) represents a physically measurable quantity, as has 
been demonstrated in several recent experiments [S71l38j . 



B. Spin-boson model and master equation 

The spin-boson model describes a two-level system 
which is linearly coupled to an environment of harmonic 
oscillators. The Hamiltonian of the model consists of a 
system part Hg, an environmental part He, and an in- 
teraction part Hj: 



H 



Hs + He 

UJQ 



Hr 



(4) 



E 



2m„ 



m n UJ n 2 



where a x and a z are Pauli matrices, loq denotes the en- 
ergy splitting of the two-level system (h is set equal to 
1), and m„, w„, x n , and p n represent the masses, fre- 
quencies, and position and momentum operators of the 
environmental oscillators, respectively. Finally, K n de- 
scribes the coupling of the system to the n-th oscillator. 

We assume that the system-environment coupling is 
weak and employ the second-order time-convolutionless 
master equation in the interaction picture to describe the 
dynamics of the open system 



dt 



P(t) = 



-i[H s (t),p(t)] (5) 



i,J=0,l 

We note that the derivation of this master equation is 
carried out without using the Markov or the rotating 
wave approximation. We have introduced the operators 
o"i = er + and do = cr_. The time-dependent coefficients 
ctij(t) form a Hermitian matrix and can be expressed in 
terms of the correlation functions of the environmental 
operator, 



m = E 



(a n e-*^ t + aU 1Unt ) 



(6) 



where a) n and a n are the creation and annihilation oper- 
ators of the environmental modes. Explicitly, we find 



aii(i) = \ I ds Re \Ti E {B(t)B(s)p E }e-^ t -^ 
* Jo L 



aoo(t) 



a w (t) 



1 



ds Re 



ds Re 



Tr E {B(t)B(s)p E }e 



Tv E {B(t)B(s)p E } 



+iu (t—s) 



-iujo(t—s) 



, (7) 
■ (8) 
(9) 
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The time-dependent system Hamiltonian Hs(t) is a di- 
agonal matrix with diagonal elements 



hi(t) = 
h {t) = 



ds Im 



ds Im 



Tr E {B(t)B(s)p E }e +l ^ t ~ s '> 
TT E {B(t)B{ S )p E }e- l ^ t - s ^ 



(10) 
(11) 



The environment is assumed to be in a thermal equilib- 
rium state initially, i.e., we have p E = Z~ 1 e~@ HE , where 
Z is the partition function and /3 = A the inverse tem- 
perature (the Boltzmann constant is set equal to 1). 

It is convenient to write the master equation ([5]) in 
terms of the Bloch vector v(t) = Trs{erp(i)}. Trans- 
forming back to the Schrodinger picture one finds pQ 



~v(t)=M(t)v(t)+b(t), 



where 



-ojq 

M(t) = | uj + a yx {t) a zz {t) 
o„(t) 



(12) 



(13) 



and b(t) = (0, 0, b z (t)^ T , with real coefficients given by 

1 /"' 

a yx (t) = -2Im{a 10 (t)} = - / ds D 1 (s) smu) s, (14) 

* Jo 



ds Di(s) coswqs, (15) 



b z (t) = an(t) - a 00 (t) = -- / dsD(s) sin w s. (16) 

* Jo 

Here, -Di(s) and -D(s) denote the noise and dissipation 
kernel of the model which can be expressed in terms of 
the spectral density J(ui, f2) of the environmental modes 
as 

D 1 (s) = 2 du J(w,fi)coth ( — J cosws, (17) 
Jo \2T J 

(18) 



D(s) = 2 / J(w, f2) sin ws. 
«/ o 

In the following we choose an Ohmic spectral density 
with Lorcntz-Drude cutoff, 



7 u Vt 2 



7T LdQ ft 2 + UJ 2 



(19) 



where f2 is the cutoff frequency and 7 defines the coupling 
strength. 



C. Numerical simulation 

We have carried out numerical simulations to deter- 
mine the non-Markovianity measure ([3} for the spm- 
boson model. For any given pair of initial states pi^(0), 
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FIG. 1. (Color online) The trace distance as a function of time 
for various temperatures (top) and cutoffs (bottom). The pa- 
rameter values of the middle curve in both plots are T — uiq, 
Q = w , 7 = O.Iojo, and vi, 2 {0) = ±(1,0, 0) T . The other 
curves show the effect of a variation of one of these parame- 
ters. 



represented by the corresponding Bloch vectors 1/1,2(0), 
one first solves the Bloch equations ( 12 1 to determine the 



behavior of the trace distance by means of 

V{t) = V( Pl (t),p 2 (t)) = 1 \^(t) - v 2 (t)\ , 



(20) 



where \vi — v 2 \ denotes the Euclidean distance. An ex- 
ample is depicted in Fig. [T] 

To obtain the total backflow of information given by 
the time integral in Eq. ([3]), we subdivide the time axis 
into sufficiently small intervals [ij,tj_|_i] and evaluate the 
quantity 



N ' = 5 E [1^+1) 



- V 2 (t i+1 )\ - - v 2 {U)\ 



where the sum is extended over those intervals in which 
the trace distance increases. The measure ([3| is then 
found by maximizing A/ 7 over all pairs of initial states. 
We have realized this maximization procedure numeri- 
cally by a Monte Carlo sampling, drawing independent 
random initial state pairs from a uniform distribution 
over the state space. 

Our goal is to study the behavior of the non- 
Markovianity as a function of parameters of the envi- 
ronment, namely the temperature T and the cutoff frc- 
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quency f2 in the spectral density. To this end, we nu- 
merically estimate the measure J\f(&) according to the 
above procedure for different values of these parameters. 
The result is a full characterization of the non-Markovian 
behavior in the spin-boson model depending on the prop- 
erties of its environment depicted in Fig. [2] For the 
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FIG. 3. (Color online) Non-Markovianity A/"($) as a func- 
tion of and T up to fi,T = IOcjo- For this range, the 
maximization has been refined with additional initial state 
pairs to yield a higher resolution. The white curve shows the 
parametrization of the resonance effect given by Eq. (231. 



FIG. 2. (Color online) Non-Markovianity Af(&) as a function 
of n and T up to Q,T = 50cjo, indicated by the logarithmic 
colorbar. 



preparation of this figure, the non-Markovianity mea- 
sure has been evaluated on a grid consisting of roughly 
5000 parameter combinations, with a denser sampling in 
the range of T, £1 < IOojo shown in Fig. [3j For each 
of these points, the Monte Carlo simulation to approxi- 
mate the maximization over all state pairs has been done 
with 12000 random initial state pairs. Comparisons and 
tests with special fixed initial state pairs suggest that the 
speckled structure in Fig. [2] is due to the maximization 
procedure and does not mirror a physical effect. We will 
now discuss and interpret these results in detail. 



III. DISCUSSION AND THEORETICAL 
INTERPRETATION 

A. Resonance effect 

Figure [2] shows the non-Markovianity defined by Eq. 
|3| as a function of temperature T and environmental 
cutoff frequency f2. Note that the values of the measure 
J\f($) are represented by a logarithmic colorbar for bet- 
ter visualization. Based on this plot, one can decide for 
which parameters the system behaves Markovian and for 
which it exhibits strong non-Markovian effects. We ob- 
serve that, roughly speaking, for small cutoffs < u>o and 
for small temperatures T < ujq the dynamics is strongly 
non-Markovian, as can be seen from the high values of 
the measure Af($>) in these regions. By contrast, for 



large values of both and T, the measure is zero and 
the dynamics is Markovian, in full agreement with the 
well-known quantum optical limit pQ. 

Our simulation results show that the non-Markovianity 
measure A/"($) features a characteristic structure: It ex- 
hibits a distinct minimum, i.e., a Markovian area within a 
non-Markovian regime for certain values of temperature 
and cutoff. This can be seen more clearly in Fig. [3] which 
shows a magnification of the range fi, T < IOujq obtained 
from a numerical simulation with higher resolution. In 
the following we develop a physical interpretation for this 
feature which obviously depends on the temperature and 
the cutoff frequency so that it cannot be explained by 
looking at the spectral density J(w, f2) only. The posi- 
tion of the minimum can be extracted from the Monte 
Carlo simulation data and the dependence is found to be 
nonlinear, contrary to first the impression. 

The environmental modes enter the master equation 
via the correlation functions in Eqs. (|7j), <[8j), and ([9]). In 
particular, the coefficients a yx (t) and a zz (t) in the master 
equation depend on the noise kernel -Di(s) given in Eq. 
( 17 1, which is determined by the effective spectral density 



J e ff(w,f2,T) 



J(a;,0)coth(^), 



(21) 



where J(u>, O) is defined by Eq. (19 1. The spectral dis- 



tribution J e ff depends on three parameters: the environ- 
mental mode frequency w, the cutoff frequency f2, and 
the temperature T. One can easily see that J e ff, regarded 
as a function of u, has exactly one maximum at some 
frequency w max and the dominant environmental modes 
are located around this frequency value. Under the res- 
onance condition uj — w max the open system sees an 
environmental power spectrum, which is approximately 
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constant around its transition frequency cuq, as is illus- 
trated in Fig. |4j In the neighborhood of the transition 
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FIG. 4. (Color online) The effective spectral density J e ff(w) 
for a fixed cutoff Q, = 1.55cjo and three different values of tem- 
perature T. The middle curve (T = 0.4tJo) shows the spectral 
density obtained from the resonance condition ( 23 1 ; the other 
two curves correspond to detuned temperature values. The 
vertical line indicates the system's transition frequency ujq. 
The bottom plot shows a magnification of the top plot. 



frequency of the system the spectral distribution is, thus, 
similar to the spectrum of white noise, and we, therefore, 
expect to find the region of predominantly Markovian be- 
havior from the condition 



d 



0. 



(22) 



This equation defines a curve Q = fl res (T) in the (Q,T) 
plane representing the points at which the system's tran- 
sition frequency ujo is exactly in resonance with the max- 
imum of the effective spectral density. Hence, the re- 
gion of Markovian behavior should be located around 
this curve. To test this prediction we use Eq. (22) to 
obtain the following resonance condition: 



^res(T) — Wo-) 



ITsmh(f) 



UJQ 



Tsinh 



v T , 



LU 



(23) 



This function is shown as a white curve in Fig. [3] By 
comparison with the dark region in the figure, represent- 
ing the region of very low or even zero non-Markovianity 



measure Af(§), we see that the resonance condition in- 
deed explains very well the emergence of Markovian be- 
havior embedded in a region of large non-Markovianity. 
We note that this picture corresponds to the results of 
Ref. [7], where the damped Jaynes-Cummings model, 
describing the interaction of a two-level system with a 
damped cavity mode, has been studied. For this model 
the non-Markovianity measure, considered as a function 
of the detuning between the system transition frequency 
and the frequency of the cavity mode, shows an analo- 
gous behavior, expressing Markovian behavior for suffi- 
ciently small detuning and a transition to non-Markovian 
dynamics above a certain threshold. 



B. Stationary rate approximation 



The master equation ( 12 ) can be derived by employ- 



ing only the Born approximation, which presupposes a 
weak system-environment coupling. The conventional 
Markov approximation consists of neglecting the time- 
dependence of the coefficients of the master equation 
(iyx(t), a Z z(t)i an d b z (t) and taking their asymptotic val- 
ues for t — > oo (stationary rate approximation). The 
resulting equation is also known as the Redfield master 
equation |48j . It provides an equation of motion with 
a time-independent generator which is not in Lindblad 
form and, hence, does not yield a completely positive 
quantum dynamical semigroup. However, for many phys- 
ical problems the Redfield equation has been shown to 
provide a useful and accurate approximation of the re- 
duced open system dynamics. 
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FIG. 5. (Color online) Non-Markovianity N(J&) as a function 
of Q. and T up to Q, T = IOojo for the master equation with 
stationary rates (Redfield equation). 



We have carried out numerical simulations to deter- 
mine the non-Markovianity measure Af(&) under the sta- 
tionary rate approximation leading to the Redfield equa- 
tion. Results are shown in Fig. [5j Quite remarkably, 



() 



the dynamics generated by the Redfield equation is still 
strongly non-Markovian in a large parameter region. Ac- 
tually, comparing Figs. [5] and [3] we observe that the over- 
all behavior of the non-Markovianity measure is quali- 
tatively and even quantitatively very similar to the case 
without stationary rate approximation. Thus, we con- 
clude that the stationary rate approximation, which is 
conventionally regarded as the Markov approximation, 
does not lead to a quantum Markov process in the sense 
of the definition for non-Markovianity used here. 



IV. CONCLUSIONS 

In this article we have quantified non-Markovian effects 
in the dynamics of the spin-boson model depending on 
properties of the environment. We have characterized the 
non-Markovianity of the dynamical map as a function of 
the temperature and the cutoff frequency in the spectral 
density. This result serves to classify different regimes in 
the parameter range. On the one hand, we can under- 
stand why the predictions of a Lindblad master equation 
fail in some applications, and, on the other hand, our 
findings may be used as a tool to devise experimental se- 
tups that feature especially large non-Markovian effects. 

During the maximization over all initial state pairs 
for the evaluation of the non-Markovianity measure, we 
found that the main contribution to the measure stems 
from the time evolution of the off-diagonal elements of 
the system density matrix, i.e., the coherences. This fact 
suggests that at least in our case, the non-Markovianity is 
mainly a quantum effect. Consequently, a description by 
means of a classical rate equation might fail to reproduce 
important memory effects. 

It is often believed that non-Markovianity is to a large 
extent only a transient phenomenon. However, our com- 
parison of the results for the non-Markovianity mea- 
sure applied to the full master equation with those for 
the master equation with stationary rate approxima- 
tion (Redfield master equation) reveals that this intu- 
ition might be misleading: Although the initial oscil- 
lations of the coefficients in the master equation decay 
over times which are short compared to the system's 
relaxation time, the non-Markovianity measure steadily 
builds up over the whole relaxation process, showing no 
special significance of the transient time evolution. Con- 
sequently, this approximation yields qualitatively similar 



results for the measure, despite the fact that it neglects 
any transient effects. 

In this context, the meaning of the standard Markov 
approximation must be revisited. The approximation of 
stationary rates in the master equation is often regarded 
as Markovian description and the rates are called Marko- 
vian rates. However, our results demonstrate that an 
approximation of this kind does not necessarily lead to 
a Markovian dynamics in the sense of our definition in 
terms of the information flow. In fact, we found that 
the master equation with stationary rates is still capable 
of describing the essential features of the non-Markovian 
behavior. It should also be noticed in this context that 
in the derivation of a Lindblad master equation the ro- 
tating wave approximation is a crucial step in order to 
obtain a completely positive quantum dynamical semi- 
group. Thus, it is the combination of rotating wave and 
standard Markov approximation that has the desired ef- 
fect of yielding a Markovian semigroup; the term Markov 
approximation can in this sense be deceptive. 



The dynamics given by the master equation ( 12 1 with- 



out the rotating wave and stationary rate approxima- 
tions not only shows strong non-Markovian behavior for 
small temperatures and/or environmental cutoff frequen- 
cies, but also features an interesting resonance effect: For 
a certain narrow range of parameters within the non- 
Markovian regime, the system exhibits Markovian dy- 
namics. We were able to explain this effect by showing 
that this Markovian behavior occurs when the system 
transition frequency is in (or close to) resonance with the 
dominant modes of the environment. This can be under- 
stood by recognizing that in this case the distribution of 
modes in the vicinity of the system frequency resembles 
a uniform distribution (white noise), which clearly favors 
Markovian behavior. The next step towards a deeper 
understanding of the resonance effect is to study the be- 
havior of the non-Markovianity measure for other mod- 
els involving, in particular, structured reservoir spectral 
densities with several maxima or spectral gaps. 
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